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1. Summary and discussion 

The density and distribution functions of Student's distribution with p degrees 
of freedom (SD P ) are given, respectively, by the formulas 

r(P±l) / T 2\-(P+ 1 )/ 2 
f P {x):= ^ 1 + - and (1.1) 

F p (x) := f f P (u)du (1.2) 

J — oo 

for all real x. Most often, the values of the parameter p are assumed to be 
positive integers. However, formula (1.1) defines a probability density function 
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for all real p > 0, and, as we shall sec, it may be advantageous, at least as far as 
proofs are concerned, to let p take on all positive real values. Let us also extend 
definitions (1.1) and (1.2) by continuity to p = oo, so that 

/oo =: if and =: $ 

are the density and distribution functions of the standard normal distribution 
(SND). 

The standard normal and Student distributions are clearly among the most 
common distributions in statistics. It is a textbook fact that the SD p is close to 
the SND when p is large, say in the sense that f p (x) — > foo(x) for each real 

p— y oo 

x. By Scheffe's theorem [8], this implies the convergence of the total variation 
distance 

1 f°° 

drv{p) = ^ J \f P ( x ) ~ f(x)\dx (1.3) 

to as p — > oo. In fact, the convergence of the SD p to the SND is presented in 
[8] as the motivating case. 

Consider also the Kolmogorov distance 

d Ko {p) ■= sup \F p (x) - $(x)| 

between the SD p and SND. It is clear that, for any two probability distribu- 
tions, the Kolmogorov distance between them is no greater than twice the total 
variation distance, and hence the convergence of the latter distance to implies 
that of the former. 

However, in the present case one can say more. For any p and q in the inter- 
val (0, oo], let g?ko(p, q) and g?tv(??7 q) denote, respectively, the Kolmogorov dis- 
tance and the total variation distance between SD p and SD g , so that d,Ko{p) = 
dKo(p,°o) and d TV (p) = d TV (p, oo). 

Proposition 1.1. 

(i) For all p and q such that < p < q ^ oo 

hdTv(p,q) = d Ko (p,q) = max (F q (x)-F p (x)). (1.4) 

2£(0,oo) 

(ii) Moreover, for each p € (0, oo) the distance d^ {p,q) is strictly increasing 
in q G [p, oo], and for each q € (0, oo] the distance c?Ko(p, q) is strictly 
decreasing in p G (0,q\- In particular, 

dKoip, q) < d Ko (p, oo) = d Ko (p) (1.5) 

for all p and q such that < p q < oo, and g?ko(p) is strictly decreasing 
in p G (0, oo]. 

Statement (ii) holds as well with d^v in place of dxo- 
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This proposition and the other results stated in this section will be proved in 
Section 2. 

The Kolmogorov distance and the total variation one are apparently the two 
most commonly used distances between probability distributions. Therefore, it 
seems natural to consider the rate of convergence of cIko{p) and, equivalently, 
<1tv(p) to as p — > oo, which is part of what is done in this paper. Actually, the 
motivation for this study comes from the discussion in [4]. In turn, the paper 
[4] was motivated by developments of [6]. 

Theorem 1.2. For any real p ^ 4 

±d TV ( P )=d Ko (p)<C/p, (1.6) 
where 



Moreover, 



C:= 1 J 7 -±^= 0.158.... (1.7) 
4 V Tre^^ 



lim pd Ko (p) = C, (1.8) 

p— >OD 



so that the constant C is the best possible one in (1.6). 
In what follows, it is assumed by default that 

a := 1/p. 

Theorem 1.2 is based on 
Theorem 1.3. For any real p ^ §| 

\d TY ip) = d Ko (p) < B(a, x a ), (1.9) 

where 

B(a, x) := (8x [5aV + a(3x 6 - 7a; 4 - 5a; 2 - 3) + 2A{x 2 + 1)] ip(x) 
+ 33a 2 (2$(a;) - 1)' 



and x a is, for any a € (0,1), the unique real root x > of the polynomial 
equation 

P(a, x) := -96(x 4 - 2a; 2 - 1) - 4a(3a; 8 - 28a; 6 + 30a; 4 + 12a; 2 + 3) 

-a 2 (20a- 4 -60a; 2 -33) = 0. (1.10) 

In fact, it will be shown (Lemma 2.6) that 

B(a,x a )<C/p (1.11) 

for all p 4. 
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Note that, since the polynomial equation (1.10) is of degree 4 in x 2 , the root 
x a can be expressed in radicals of polynomials in a. 

By the triangle inequality, (1.6) implies that i d,Tv(p,q) = dKo(p,q) ^ 
^Ko(p) + ^Ko(?) < C/p + C/q for any real p and q that are no less than 4. 
Taking (1.4) and (1.5) into account, one sees that (1.9) and (1.11) immediately 
yield better bounds: 



Corollary 1.4. For all p and q such that 4 ^ p < q ^ oo 

hdrvip, q) = d Ko (p,q) < B(a,x a ) < C/p. 



(1.12) 



Graphs of the bounds B(a, x a ) and C jp are shown in Figure 1, along with the 
corresponding graph of d^ Q {p). This is done for the values of p £ [||, 30], even 
though the upper bound C/p on c?Ko(p) has been established only for p 4. 




10 



15 



P 



20 25 30 

Fig 1. Bounds B(a,x a ) (red) and C/p (blue), compared with d^o(p) (black). 



The relative errors "^'(p/ _ ^ an< ^ <* «w — 1 of the bounds in (1.9) and (1.6) 
are shown in Figure 2, for p £ [1,3.95] in the leftmost panel, for p £ [3.95,4.05] 
in the middle panel, and for p £ [4.05, 30] in the rightmost one. 
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Fig 2. Relative errors ~^~t~t — 1 (red) and d — 1 (blue) of the bounds in (1.9) and 
(1.6). 
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It appears that the bound C/p would be more accurate than B(a,x a ) for 
p G [1,3.98]; remember, however, that the bound B(a,x a ) was established only 
for p ^ || , from which the bound C/p was deduced only for p ^ 4. Anyway, the 
smaller values of p > may be of lesser interest, since for such p the Student 
distribution is not very close to the standard normal one. On the other hand, 
for large p the bound B{a 1 x a ) appears significantly more accurate (in terms of 
the relative errors) — but much more complicated — than the bound C/p. Yet, 
even for p as small as 4, the relative errors of the bounds C/p and B(a, x a ) are 
both only about 5%, with the corresponding absolute errors less than 2 x 10~ 3 . 
For p = 12, the relative and absolute errors of the bound C/p are less than 1.5% 
and 2 x 10~ 4 , respectively, and the corresponding figures for the bound B(a, x a ) 
are about 0.5% and 6 x 10 -5 . Also, by (1.8), the relative error d ^ ^ — 1 of the 
upper bound C/p goes to as p — > oo; in view of (1-11), the same holds for 
the upper bound B(a, x a ). One may as well note that, if the distance cZko(p) is 
considered as a kind of "initial" error — of the approximation of the Student 
distribution by the SND, then the relative error , . — 1 is a relative error "of 
the second order" , in the sense that it is the relative error of the estimate C /p 
of the initial error cZko(p); the same statement holds with B(a,x a ) in place of 
C/p. 

Figure 2 also suggests that the threshold value 4 in the condition p ^ 4 in 
Theorem 1.2 is very close to the best possible one for which the comparison 
(1.11) between the bounds in (1.6) and (1.9) is still valid. 

In [7], an asymptotic expansion for the tail 1 — F p (x) of the SD P was obtained, 
which provides successive approximations (say A p j(x)) that are good for very 
large values of x, as illustrated in Figure 3 — for p = 14. The right panel 




Fig 3. Left panel: the successive approximations A\4 y i (x), . . . ,Ai4 t a(x) of 1 — i<i4(a;) as in 
[7], colored red, orange, green, cyan, blue, purple, respectively; the graph of the tail 1 — Fi4(x) 

is black. Right panel: the graphs of the corresponding absolute relative errors \ ~ 1 

(j = l,...,6). 

of Figure 3 suggests that each approximation A p j(x) has its own "maximum 
competency" zone of values of x, for which it is the best, over all j's; it appears 
that this zone is a neighborhood of 00, which gets narrower as j increases. 
Clearly, the bounds given in the present paper differ quite significantly in kind 
and purpose from those given in [7]. 
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2. Proofs 

The main idea of the proof of the inequality in (1.6) is to reduce it, through a 
number of steps, to systems of algebraic inequalities. Such systems, by a well- 
known result of Tarski [9, 1] (rooted in Sturm's theorem), can be solved in a 
completely algorithmic manner. Similar results hold for certain other systems 
which may also involve the logarithmic function (whose derivative is algebraic), 
the SND density function ip (whose logarithm is algebraic) , and the SND distri- 
bution function $ (whose derivative is ip). The bound B(a,x a ) in Theorem 1.3 
is such an expression. The Tarski algorithm is implemented in latter versions of 
Mathcmatica via Reduce and other related commands. For instance, a command 
of the form 

Reduce [condl kk cond2 kk {varl,var2,. . . ,}, Reals] 

returns a simplified equivalent of the given system (of equations and/or in- 
equalities) condl. cond2, . . . over real variables varl, var2, .... However, the 
execution of such a command may take a very long time (and/or require too 
much computer memory) if the given system is more than a little complicated, 
as is e.g. the case with the system B(a,x a ) < C/p & a = 1/p & p ^ 4, which 
provides the way to deduce the bound in (1.6) from that in (1.9). Therefore, 
Mathematica will need some human guidance here. It appears that all such cal- 
culations done with the help of a computer are, at least, as reliable and rigorous 
as the same calculations done only by hand. 

The main difficulty to overcome in this paper was to construct the upper 
bound B(a,x a ) on c?Ko(p), which would be, on the one hand, accurate enough 
and, on the other hand, provide a traversable bridge from d,Ko(p) to the simple 
upper bound C/p, as indicated above. In turn, the bound B(a,x a ) was ob- 
tained in several steps, described in detail in the statements of Lemmas 2.1-2.5, 
presented in Subsection 2.1. 

The first step is to note that the difference foo(x) — f p (x) between the densities 
of the SND and SD p changes its sign exactly once, from + to — , as x increases 
from to oo (Lemma 2.1). A key observation here (essentially borrowed from 
[5]) is that, luckily, for the defined in (2.1) ratio r p (x) of the densities, the 
logarithmic partial derivarive lnr p (x) increases in x £ [0, 1] and decreases in 
x £ [1, oo) — with the same switch-point 1 for all p > 0. This implies that r p {x) 
decreases in x <S [0, 1] from r p (0) < 1 and then increases in x S [1, oo) to oo, so 
that the difference i 7 ^ — F p between the SND and SD p distribution functions 
switches its monotonicity pattern just once — from increase to decrease - 
on the interval [0, oo), which provides a more manageable expression for the 
Kolmogorov distance cIko(p)- 

The next step concerns the difficulty that the expression (1.1) for f p (x) con- 
tains the so-called Wallis ratio F( £ ii)/F(|), which is not algebraic, and whose 
logarithm or derivative or logarithmic derivative is not algebraic either. To deal 
with this problem, we have just developed in [3] series of high-precision upper 
and lower algebraic bounds on the Wallis ratio; for the purposes of the present 
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paper, the first upper bound and the second lower bound in the corresponding 
series in [3] already suffice (Lemma 2.2). (A recent paper [2] provided other 
new upper and lower bounds on the Wallis ratio, improving on a number of 
preceding results. The series of bounds given in [3] (except a few first members 
of those series) are tighter than all the bounds in [2].) By using the mentioned 
results of [3], we obtain an upper bound, written as H(a, x)/\/2n, on the differ- 
ence foo{x) — fp(x) between the densities of the SND and SD p , which has an 
algebraic expression in place of the Wallis ratio (Lemma 2.3). 

According to (2.3), d,Ko{p) equals a definite integral (in x) of the difference 
/oo (%) — fp (z) ; so, this integral can be bounded from above by the corresponding 
integral of the just mentioned upper bound H (a, x) / \/2ir. However, the latter 
integral is still problematic to estimate accurately enough. Toward that end, 
by some tweaking of the third-order Taylor polynomial in a for H (a, x) near 
a = 0, we construct an upper bound H2(a,x) on H(a, x), which is just the 
product of (f{x) and a polynomial in a, x (Lemma 2.4). Thus, the bound Hi{a, x) 
has certain nice properties (Lemma 2.5). Also, the relevant definite integral of 
H2{a, x) (corresponding to the mentioned one of H(a, x)) can be easily expressed 
in terms of the functions ip and $, thus finally resulting in the bound B(a,x a ) 
in (1.9). 

Inequality (1.11) (which, together with (1.9), yields the inequality in (1.6)) 
is provided by Lemma 2.6, whose proof is rather technical and relies on the 
Mathcmatica command Reduce, as described above. As for Proposition 1.1, it 
follows easily from Lemma 2.1 and the result of [5]. 

It appears that essentially the same method can be used to obtain even tighter 
upper (as well as lower) bounds on the distances cZko(p) and cZtv(p); toward 
such an end, one could use bounds in [3] on the Wallis ratio of higher orders 
of accuracy, as well as tweaked- Taylor polynomials for H (a, x) of higher orders. 
The limitations on the attainable accuracy of such bounds on cIko(p) appear to 
be mainly set by the existing computational power; also, the proofs of the yet 
tighter bounds can be expected to be even more complicated. 

In accordance with the above description of the scheme of proof, the current 
section is organized as follows. In Subsection 2.1, the mentioned lemmas are 
stated, thus presenting most of the main steps of proof. Next, in the same 
subsection, Proposition 1.1 and Theorems 1.3 and 1.2 are proved based on these 
lemmas. Finally, in Subsection 2.2 the lemmas stated in Subsection 2.1 (and 
requiring proof) arc proved. 



2.1. Statements of lemmas, and proofs of the main results 

Introduce 

r P ,q(x) := yj^- and r p (x) := r p>00 (x) = ^My- (2.1) 

Jq\ x ) l P\ x ) 

Lemma 2.1. For each pair (p, q) such that < p < q ^ oo 

(i) the ratio r Pj9 (x) decreases in x € [0, 1] from r Ptq (0) < 1, and then increases 
in x € [1, oo) to oo; therefore, 
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(ii) there is a unique point x Pjq G (0,oo) (which is in fact greater than 1) such 
that 

fp{x) < f q (x) for all x G [0,x Ptq ), 

fp{ x p,q) = fq( x p,q)> (2-2) 

f P (x) > f q (x) for all x G (x M ,oo), 



and hence 
For brevity, let 



dKo(p,q) = F q (x p , q ) - F p (x Pt g). (2.3) 
Xp • — Xp^QQ. (2.4) 



Lemma 2.2. For all real p > 

L 2 (a) < r„(0) < Utia), (2.5) 

where 

This follows by the main result in [3]; the notations r p (0), Lk(a), and Uk(a) 
in the above Lemma 2.2 correspond to r(p), Lk{p), and Uk(p) in [3]. 

The first inequality in (2.5), together with the definition (1.1), immediately 
yields 

Lemma 2.3. For all real p > and all x G K 

, / x , / x H{a,x) 

foo{x) ~ fp{x) < 



/2tt 
w/iere 

H(a,x) := er x2 / 2 - L 2 (a)(l + ax 2 )-^ . 

By some tweaking of the third-order Taylor polynomial in a for H (a, cc) near 
a = 0, one obtains 

Ma, X ):= a -^e-^, (2.6) 

where P(a,x) is as in (1.10), so that H2(a,x) be an upper bound on H(a,x): 
Lemma 2.4. For all a G (0, ||] and x G (0, ^) 

H(a,x) < H 2 (a,x). 

Lemma 2.5. For each a G (0, 1), there is a unique real root x > of the 
polynomial equation (1.10), so that i a is correctly defined in the statement of 
Theorem 1.3. Moreover, 

H2(a, x) > for all x G (0,x a ), 
H 2 {a,x a ) = 0, 
H 2 {a,x) < for all x > x a . 
Furthermore, x a is strictly and continuously increasing in a € (0, 1). 
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Lemma 2.6. For all p ^ 4 inequality (1.11) holds. 

Now one is ready to prove Proposition 1.1 and Theorems 1.3 and 1.2, which 
will be done in this order. 

Proof of Proposition 1.1. Take indeed any p and q such that < p < q ^ oo. 
By Lemma 2.1 and the symmetry of the SD P , 



d T vM = / (/,-/*)+ / ih-h) 

J 0Cp,q 

2 / P '\f q - f P ) = 2(F q (x M ) - F p {x p , q )) = 2d Ko (p,q), 





which proves part (i) of Proposition 1.1. Part (ii) of the proposition now follows 
by the second equality in (1.4) and the stochastic monotonicity result of [5], 
which implies that F p (x) is strictly increasing in p £ (0, oo] for each x £ (0, oo). 

□ 

Proof of Theorem 1.3. The equality in (1.9) immediately follows from Proposi- 
tion 1.1. Take any a £ (0, ||] (corresponding to p ^ ||). We claim that x v < x a , 
where x p and x a are as in (2.4) and Lemma 2.5, respectively. Assume the con- 
trary, that x p > x a - Note that H2(§§ ; t|t) < 0; so, by Lemma 2.5, £29/50 < i|f 
and hence x a < ^ for all a £ (0, ||]. Therefore, in view of Lemmas 2.3 and 
2.4, 

V2^(foo(x) - f P (x)) < H{a,x) < H 2 {a,x) (2.7) 

for all x £ (0, x a ] — still assuming that a £ (0, On the other hand, by 
Lemma 2.1, ^ foo{x) — fp(x) for all x £ (0, a; p ] and hence, by the assump- 
tion, for all x £ (0, x a ]. Now (2.7) implies < H2(a,x a ), which contradicts 
Lemma 2.5. Thus, indeed x p < x a . Recalling now (2.3) and using (2.7) and 
(again) Lemma 2.5, and also recalling (2.6), one has 

rx p 

d-Ko(p) = Foc(%p) - F p (x p ) = / (foo(x) - f p (x)) dx 

</ max <r m^. 



2tt Jo V2tt 

It remains to verify that du = B(a,x), which can be done either by 

hand or using Mathematica. The proof of Theorem 1.3 is now complete, modulo 
the lemmas. □ 

Proof of Theorem 1.2. The relations in (1.6) immediately follow by Theorem 1.3 
and Lemma 2.6. It remains to verify (1.8). First here, use l'Hospital's rule to 
find that for all real x 

lim = Urn %Q = X (x) := *- - ^ = 1 <p(*)i (™) 

aio a aio da 4 
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the second equality in (2.8) can be obtained either using the Mathcmatica com- 
mands D (for differentiation), Simplify, and Limit or otherwise. 
Next, introduce 



c a := /i/a(0) and g a (x) = fi/ a (x)/c a 



(2.9) 



for all real a ^ 0, assuming the convention 1/0 := oo, so that fi/ a (x) — c a g a (x). 
Then for all real a ^ and all real x 

\h/a{x)- fryo{x)\ < \c a -c \g a (x) + c \g a (x)-g (x)\ ^ \c a -c \ + \g a (x)-g (x)\, 

(2.10) 



\Cg— go] 



IS 



since g a {x) ^ 1 and cq = 1/V27T < 1. By (2.8) and (2.9), the ratio 
continuous in a > and converges to a finite limit (</j(0)/4) as a ], 0, and hence 
is bounded in a G (0, 1]. Now note that 



dg a {x) 



da 



where 



(Dg)(a,x) := 



(l + ax 2 )- {1+3am2a) \(D9)(a,x)\ «S \(Dg)(a,x)\, 
(l + ax 2 ) In (l + ax 2 ) - a(l + a).T 2 



2a 2 



Using the Taylor expansion ln(l + u) = u — 9u 2 /2 for u > and some 6* = 6*(u) G 
(0, 1), one sees that (Dg)(a,x) is a polynomial in a,x,9 and hence bounded in 
(a,x) € (0,1] x [0, xo] — note that, in accordance with the definition of x a in 
Theorem 1.3, 



x 



= A/l + ^e (0,oo); 



hence, 



da 



is bounded in (a,x) € (0,1] x [0,xq] and, by the mean value 



theorem, so is 



\ga{x)-ga{x)\ 



Recalling also (2.10) and that the ratio 



\Cg-C \ 



IS 



I f ! \ f ( \ I 

bounded in a £ (0, 1], one concludes that the ratio 1/a — °° is bounded in 
(a,x) e (0,1] x [0,xo]. So, by (2.8) and dominated convergence, 

fl/a(x) ~ foo(x) 



pd Ko (p) > P [Foo(xq) - F p (x )] 



aiO 



■ dx 
a 



where A(x) is defined in (2.8). This, together with (1.6), implies (1.8). The proof 
of Theorem 1.2 is now complete, modulo the lemmas. □ 



2.2. Proofs of the lemmas 

Proof of Lemma 2.1. Take indeed any p and q such that < p < q oo. A key 
observation here (borrowed from [5]) is that r p , q (x) decreases in x G [0, 1] and 
increases in x G [1, oo). Moreover, by the lemma in [5], / p (0) increases in p > 
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and hence r Pi9 (0) < 1. On the other hand, it is easy to see that r Ptq (x) — > oo 
as x — > oo. This completes the proof of part (i) of Lemma 2.1, which in turn 
implies that there is a unique x PA £ (0, oo) such that r Pt g(x) < 1 for x £ [0, x Piq ), 
r p,q( x p,q) = 1; an d r p,q{ x ) > 1 f° r x € oo) (at that necessarily x p>q > 1). In 

other words, one has the relations (2.2). Since (F q — F p )' = f q — f p , one now sees 
that F q (x) — F p (x) increases in x £ [0,x Pig ] from to F q (x p . q ) — F p (x p . q ) > 0, 
and then decreases in x £ [x Ptq ,oo) to 0. So, (2.3) follows by the symmetry of 
the Student and standard normal distributions. Thus, the lemma is completely 
proved. □ 

Proof of Lemma 2.4- Indeed assume that a £ (0, ||] and x £ (0, ^ ). Consider 
the difference 

5(a) := 5(a, x) := H(a, x) - H 2 (a, x) = ?Jf' x *} 2 - L 2 (o)(l + ai 2 )"^, 

384e 2: ' z 

where 

P(a,x) := 384- aP(a,x). 

We have to show that S(a, x) < 0. Obviously, the system of inequalities P(a, x) ^ 
0, < a ^ and < x < ^ is algebraic and thus, by the well-known result 
of Tarski [9] can be solved completely algorithmically. The Mathematica com- 
mand Reduce [tP<=0 && 29/50>=a>0 && 123/50>x>0] outputs False, where 
tP stands for P(a,x). This means that P(a,x) > — for all a £ (0, ||] and 
x £ (0, 4|f )■ So, 5(a, x) equals 

5(a) := 5(a, x) := In In (i 2 (a)(l + a* 2 )"^) 

in sign. Introduce 

/7-»r\/ \ a 2r// \ 4aQ(a,x) 2(1 + a) - 2 ^ 

(D5)(a) := 4a 5 (a) = - i- - 2 In (1 + ax ) 

V A ' W P(a,x) 1 + ax 2 V ' 

If 21 6 1 \ 

+ 6( v 72a + TT^-TT^ + iT^- 4 J' 

(DDS)(a) := ( ^ 3 (a) (1 + a) 2 (l + 2a) 2 (l + 3a) 2 (l + ax 2 ) 2 P(a, x) 2 , 
where 

Q(a, x) := a 3 (20x 4 - 60x 2 - 33) - 96a (x 4 - 2x 2 - l) - 768. 

Note that (DD5)(a) is a polynomial in a and of degree 11 in a and of degree 
20 in x. The command Reduce [DDde>=0 && 29/50>=a>0 && 123/50>x>0] out- 
puts False, where DDde stands for (DD5)(a). This means that (DD5)(a) < 
— for all a £ (0, and x £ (0, ^j-)- On the other hand, one can check (using 
Mathematica or otherwise) that (D5)(0+) = 5(0+) = 0. Thus, one concludes 
that indeed 5(a, x) < and hence 5(a, x) < 0. □ 
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Proof of Lemma 2.5. Take indeed any a G (0,1). By (2.6), H2(a, x) equals 
P(o, x) in sign. So, the first two sentences of Lemma 2.5 can be proved using the 
Mathematica command Reduce [P>0 && 0<a<l && x>0,x]. That x a is strictly 
increasing in a G (0,1) now follows by the command Reduce [PP [a, x] ==0 
==PP[b,y] && 0<a<b<l kk 0<y<=x] , which (takes about 15 seconds on a stan- 
dard laptop and) outputs False; here, PP[a,x] stands for P(a,x). Finally, the 
continuity of x a in a can be verified by the implicit function theorem; here, it is 
enough to check that §^(a, x a ) ^ for all a G (0, 1), which can be done using 
the command Reduce [P==0 kk DPx==0 kk 0<a<l kk x>0] (with DPx standing 
for j^-(a,x)), which outputs False. □ 

Proof of Lemma 2.6. By Lemma 2.5, a i— > x a is a one-to-one map of (0, j] onto 
(So, £1/4]. Let (xo,Xiu] 3 x i-» a x G (0, j] be the corresponding inverse map. 
So, it suffices to show that B(a x ,x) < Ca x for all x G {xo,Xiu]. Assume indeed 
that x G {xQ,Xi/i\ and consider the ratio 

p{x) := l!i "'- n :i (2.11) 

a x 

Introduce also 

qi{x) :=3 + 12a; 2 + 30a; 4 - 28a: 6 + 3a; 8 , 

q 2 (x) := - 783 - 2952a; 2 - 1284a; 4 + 2952a; 6 - 234a; 8 - 1608a; 10 

+ 964a; 12 - 168a; 14 + 9a; 16 , 
q 3 (x) :=33 + 60a; 2 - 20a; 4 . 

The command Reduce [q3 <= kk xxa [0] < x <= xxa [1/4]] (with q3 and 
xxa[a] standing for qs(x) and x(a)j outputs False, which shows that q^{x) > 0. 
Now using the command Reduce [P==0 kk 0<a<=l/4 kk xxa [0] <x<=xxa [1/4] ] , 
where P stands again for the polynomial P(a, x) as in (1.10), one finds that 



Qi(x) + \/q2{x) 
Q3(x) 



moreover, a x > and qs(x) > imply that qi(x) + y q 2 (x) > 0. So, in view of 
(2.11), 



e - 2 /2 



p'(x)—— 2^^q 2 {x)( qi (x) + y / q 2 {x)Y 

= V2(pao(x) +poi(x)\ / q 2 (x) ) + e x /2 Cy/n(pw(x) +pu(x) y/ q 2 {x) ), 
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where 

p m {x) :=105705x - 1945539a; 3 - 13305006a; 5 - 26650971a; 7 - 3174714a; 9 
+ 49512627a; 11 + 23388786a; 13 - 45078003a; 15 - 9879213a; 17 
+ 26892909a; 19 - 5379786a; 21 - 8094383a; 23 + 6008972a; 25 

- 1844301a; 27 + 296622a; 29 - 24 4 35a; 31 + 810a; 33 , 
p m {x) := - 21789a; - 259929a; 3 - 492804a; 5 + 366741a; 7 + 967263a; 9 

- 120468a; 11 - 487080a; 13 + 188214a; 15 + 177266.x 17 - 151973a; 19 
+ 43674a; 21 - 5625a; 23 + 2 70a; 25 , 

p w (x) := - 1368576 + 9287136a; 2 + 67830048a; 4 + 113324832a; 6 - 52129440a; 8 

- 230541408a; 10 + 74263392a; 12 + 151161696a; 14 - 110996640a; 16 
+ 30085440a; 18 - 3715200a; 20 + 172 8 00a; 22 , 

Pll {x) :=171072 + 1974240a; 2 + 2340576a; 4 - 4409568a; 6 - 3045600a; 8 
+ 2911680a; 10 - 700800a; 12 + 57600a; 14 . 

Executing now the command Reduce [z<=0 && xxa[0] <x<=xxa[l/4] ] with z 
standing for pio(x) + pn(x)y / q 2 (x), one sees that pio(x) + pn(x)y / q 2 (x) > 0, 
so that p'(x) equals 

V2(p m (x) + p i(x)^J q 2 {x)) 

Pl(x) := ~^JTi — TTT — , \ + c ^ 

e x I 1 Xpva(x) + pn(x)^q 2 (x)) 



m sign. 
Next, 



P2 (x) :=p[(x)e x /2 {pio{x) + Pll {x)^q 2 (x)) 2 

=Cq(x) + ci(x)yj q 2 (x) + c 2 (x)q 2 (x) + c 3 (x) / ^ q 2 (x) , 



where 



c (x) :=p 10 (x)p' Q0 (x) - p 00 (x)p' w (x) ~ xp 00 (x)p 10 (x), 

d(x) :=p n (x)p' 00 (x) - p 00 (x)p' n (x) - xp 00 (x)p n (x) + p 10 (x)p 01 (x) 

- Poi(x)p[ (x) - xp i(x)pio(x), 
c 2 (x) :=p u (a;)poi(a;) - Poi(x)p' n (x) - xp i(x)p n (x), 
c 3 (x) :=q 2 (x)(p 01 (x)p 10 (x) ~ p 00 (x)p 11 (x)) /2. 

The command Reduce [rho2 [x] <=0 kk xxa [0] <x<=xxa [1/4] ] outputs False. 
So, p 2 (x) > (for all x <E (xo,x 1 / 4 \) and hence pi(x) increases in such x. 
Moreover, pi(xo) = 0, which implies that p\ > 0. That is, p' > and p is 
increasing on the interval (icii/J, to p(x\/4) < 0. Thus, p < on (Jo, £1/4], 
which implies that indeed B(a x ,x) < Ca x for all x £ (xo,Xi/i\. □ 
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